arXiv:hep-ph/0211299v3 16 Feb 2003 


TU-675 
Nov. 2002 


Top mass determination and 0(ct 5 s m) correction to 

toponium IS energy level 


Y. Kiyo and Y. Sumino 


Department of Physics, Tohoku University 
Sendai, 980-8578 Japan 


Abstract 

Recently the full 0(a 5 s m, log ag) correction to the heavy quarkonium 15 energy 
level has been computed (except the a 3 -term in the QCD potential). We point out that 
the full correction (including the log ag-term) is approximated well by the large-/?o ap¬ 
proximation. Based on the assumption that this feature holds up to higher orders, we 
discuss why the top quark pole mass cannot be determined to better than O(Aqcd) ac¬ 
curacy at a future e + e~ collider, while the MS mass can be determined to about 40 MeV 
accuracy (provided the 4-loop MS-pole mass relation will be computed in due time). 
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Recently a large part of the 0(a 5 s m) corrections [IJ ;2j to the energy spectrum of the 
heavy quarkonium IS 1 state has been calculated. Combining this with the previously known 
0(a F mlogas) corrections [310] the only remaining piece to be computed in order to complete 
the 0{a%m) corrections is the non-logarithmic term ( 03 ) of the static QCD potential at 3-loop. 
Using a Pade estimate JSJ of 03 , Ref. [2] examined the scale dependences and the convergence 
properties of the bottomonium 15 and the (would-be) toponium 15 energy levels. The depen¬ 
dences of the energy levels on the value of a 3 are found to be rather weak. As for the toponium 
case, Ref. j2j concluded that the top quark pole mass can be extracted from the 15 energy 
level with a theoretical error of about 80 MeV. This estimate of the theoretical error on the 
top quark pole mass appears to be considerably smaller as compared to a previous common 
consensus that the pole mass has a theoretical uncertainty of order Aq CD ~ 200-300 MeV [ 6 j. 

In this paper we discuss two issues. First we point out that the presently known 0{oP s m) 
correction to the 15 energy level is approximated fairly well by its large-/?o approximation (naive 
nonabelianization) [7j. We consider this fact to be quite non-trivial because of the following 
reason. We know that from 0(a 5 s m ) the ultrasoft scale starts to contribute to the energy level. 
Since it is a completely new type of contribution (as compared to the lower-order corrections), 
and since it is generally believed to give very large corrections M, we have expected that 
the large-/3o approximation may well fail to be a good approximation at 0(a 5 s m ) in the energy 
level. 

One may wonder that our point, that the large-/?o approximation is good, is in contradiction 
to the conclusion of [SJ: “We have found that the N 3 LO corrections are dominated neither 
by logarithmically enhanced a|ln(a,s) nor by the renormalon induced /3 q«| terms and thus 
the full calculation of the correction is crucial for quantitative analysis.” I 11 fact, there is no 
contradiction, because the definition of the “/3o a i terms” in JZJ differs from that of the usual 
large-/?o approximation.* Nevertheless, we have to say that the above statement of J2, is quite 
misleading, since it does not address the difference between its terms and the large-/3o 

approximation, and since it is the large-/3o approximation that is the empirically successful 
approximation and, therefore, the renormalon dominance picture has often been discussed in 
this context in the literature. 

Secondly, we discuss an error estimate of the top quark pole mass based on the assumption 
that the large-/3o approximation continues to be a good approximation up to higher orders. At 
the same time we discuss the accuracy with which the top quark MS mass can be extracted 
from the toponium 15 energy level. 

The sum of the full 0(a 5 s m ) and 0(oigmlogas) corrections to the energy level of the heavy 
quarkonium 15 state is given in Eqs. ( 6 ), (12) and (13) of (2j. The part unrelated to the lower- 
order corrections via the renormalization-group equation (for the running of the coupling) can 
be extracted by setting L /x = log [///(Giro's (/i)m po i e )] = 0. It reads numerically 


5E[ 3) 


{C F ois(n))‘ 


( <xs(v) \ 


c 3 ~ 7078.8 + 0.03125 a 3 


^pole X ^ J C 3 , 

1215.5 n f + 69.451 n f 2 


1.2147 n f 3 + 474.29 log (a s (/x)) , 


( 1 ) 

( 2 ) 


*For instance, the term proportional to cl\0q is not included in the /3q term of |2], whereas a part of ai/?o is 
included in the large-/3o approximation. 
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a s 

0.1 

0.2 

0.3 

C3 | rif =4 

2354 

2683 

2875 

t-3 |rif=5 

1613 

1942 

2134 

t-3 | n,y=4, log(ag)—>0 

3446 

3446 

3446 

C 3 |nf=5, log(as)—>0 

2705 

2705 

2705 

c 3 (large-/3 0 ) \ nf =4 

2456 (104%) 

2456 (92%) 

2456 (85%) 

c 3 (large-#) \ nf=5 

1913 (119%) 

1913 (98%) 

1913 (90%) 


Tabic 1: Numerical values of C 3 and its large-/?o results are shown for as = 0.1, 0 . 2 , 0.3 and rif = 4, 5. For 
each C 3 (large-/ 3 o), the ratio to the full result (C3) is shown in the parenthesis. The Pade estimate 0 of 03 is 
used in Eq. m to obtain C3. 


where Cp = 4/3 is a color factor. 

In general, the large-#) approximation of a quantity, at a given order of perturbative ex¬ 
pansion in as, is defined as follows: We first compute the leading order contribution in an 
expansion in 1 /n/, where rif is the number of light quark flavors, which comes from so-called 
bubble chain diagrams. Then we transform this large rif result by a simplistic replacement 
n f n f ~ 33/2 = —(3/2)#j. I 11 many phenomenological applications the large-#) approxi¬ 
mation turns out to be a good approximation of the full result for quantities which contain 
the leading renormalon, see e.g. [El El CHI HU- The corresponding correction to Eq. m in the 
large-#) approximation is given by dim 


c 3 (large-#)) 




3/517 19 7T 2 11C 3 7T 4 7T 2 C 3 3Cs\ 
Po \864 144 6 1440 8 2 ) 

5649.36 - 1027.16 n f ± 62.2519 n f 2 - 1.25761 n f 3 


(3) 


In Table [H we compare C 3 and C 3 (large-#)) for values of rif and as corresponding to the 
Y(IS') and toponium IS' states. For a 3 , we used the Pade estimate H as well as the estimate 
based on the renormalon dominance picture fH]; c 3 differs by less than 3% when we use these 
estimates, for rif = 4, 5/ We also varied a 3 by ±100% in Eq. (J2J) and find that c 3 changes by 
less than ± 10 % for rif = 4, 5. As we can see from the table, the large-#) approximation turns 
out to lie between 85% and 120% of the full result in the relevant cases. We observe that the 
agreement becomes substantially worse if we remove the log as term from the full result. 

In Figs. Ct) and b), we show the renormalization scale (/i) dependences of the IS energy 
level when we use the pole mass and the MS mass*, respectively, to express the energy level. We 
used the e-expansion P3 to cancel renormalons in the MS mass scheme; the relevant formulas 
are given in the Appendix. Fig. GJi) is essentially a reproduction of Fig. 2(b) of [2], by including 
the leading order (LO) curve in addition. As pointed out by |2], the next-to-next-to-next-to- 
leading order (NNNLO) prediction becomes insensitive to the scale variation at p ~ 15 GeV, 

Mhe corresponding estimates of the three loop coefficient 03 are given by as(Pade)/4 3 = 98, 60, and 
a 3 (Pineda)/4 3 = 72, 37 for rif = 4, 5, respectively jSJ 77- 

Mhe pole-MS mass relation is known up to 3 loops presently. The 4-loop correction is replaced by its large-/?o 
approximation in our analysis. 
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Figure 1: The renormalization scale dependences of the energy level of the vector toponium 15 state for a) 
the pole mass scheme and b) MS mass scheme, respectively. The solid curves are the NNNLO results, dotted, 
dashed and dot-dashed curves denote the LO, NLO and NNLO results, respectively. 


and that the sum of the 0{oP s m) and O(a^mlogas) corrections becomes small around this 
scale. On the other hand, in Fig. nr>), we see a good convergence behavior at fi ~ 50-80 GeV. 
In Fig. [21 the vertical scale is magnified and the scale dependences of the energy levels at the 
NNNLO in both mass schemes are compared. We find a much better stability of the prediction 
in the MS mass scheme over a wide region 40 GeV < /jl < 160 GeV. From this analysis, we 
consider the scale /x ~ 50-80 GeV to be an optimal scale choice in the MS mass scheme. By 
varying // between 30-160 GeV, we estimate the theoretical error of the MS mass to be order 
40 MeV at NNNLO if it is extracted from the IS energy level. (We obtain an error of about 
200 MeV if a similar estimate is applied for the pole mass.) 

In the pole mass scheme, it is natural to choose the renormalization scale around the Bohr 
scale, /x ~ Cfc*s m ~ 30 GeV. This is because there is only one logarithm \og[ii / {CfO-stti)] in 
the energy level, associated with the renormalization scale /xp and because this logarithm is 
minimized around the Bohr scale. On the other hand, in the MS mass scheme, two types of 
logarithms log (/i/m) and log ^/(F^asm)] are included in the expression for the energy level,^ 
where m = "^ms( m Ms) the renormalization-group invariant MS mass; see Eqs. (Il3|) , (1171) and 
m■ Therefore, a natural scale, which minimizes the logarithmic contributions, lies between 
the Bohr scale and the hard scale, CfOisiti < /x < m. This aspect of the renormalization 
scale, when the leading renormalon uncertainty is removed, has been discussed already for the 
bottomonium energy levels HE and a further detailed study of the scale choice (in the context 
of the QCD potential) has been given in [ TTj . 

If we replace c 3 by c 3 (large-/3 0 ), the corresponding figures to Figs. Hk,b) and El look very 

§At NNNLO, the logarithm associated with the ultrasoft scale ~ a 2 s m is not accompanied by the renormal¬ 
ization scale /x. 

^This stems from the fact that one needs to expand the pole mass and the binding energy in the same 
coupling constant as(/x) in order to achieve the decoupling of infrared degrees of freedom at each order of the 
perturbative expansion. 
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Figure 2: The IS 1 energy levels of the vector toponium state at NNNLO are plotted in the pole and MS 
schemes. A horizontal line, Mis = 347.4 GeV, is drawn for a guide. 

similar; these were shown in UB[. The main observations in the analysis in the large-/?o approx¬ 
imation were as follows (1) hi the pole mass scheme, with any choice of the scale p, 

the perturbative series of the IS 1 energy level does not show a healthy convergence behavior, 
hence the level cannot be predicted with an accuracy better than O(Aqcd)- (2) In the MS mass 
scheme, one observes a good convergence of the perturbative series in the range m« s </i<m, 
as well as stability of the prediction in this range. Both of these observations still hold at the 
best of our present knowledge. It is intriguing whether these features will remain valid even 
when a 3 and the 4-loop relation between the pole and MS masses are computed fully in the 
future. 

Let us address how the theoretical error of about 80 MeV for the top quark pole mass was 
obtained in Ref. |2J. It is dominated by the uncertainty induced by the error of the input 
as(M z ). The uncertainty due to the scale dependence was estimated by varying // between 
10-30 GeV and an error of 20.5(=41/2) MeV was assigned as an uncertainty from this source. 
Uncertainties from other sources were estimated to be even smaller. Here, let us concentrate 
on the error estimate from the scale dependence and discuss its problem. The smallness of 
this error ensures, partly, the smallness of the total error (80 MeV). However, if the same 
estimation method is applied to the LO and NLO curves in Fig. Ek), we should infer that 
the NNLO correction is small, in contradiction to its true large size. Thus, apparently there 
is a danger in relying on this estimation method. By contrast, our error estimate of the top 
quark MS mass from the scale dependence does not suffer from the same problem. The same 
estimation method works at lower orders, because the perturbative series in Fig. shows a 
healthy convergence behavior and the scale dependence decreases as we include more terms 
around the relevant scales. 

At this stage, there seems to be a puzzling point: on the one hand, the validity of the large- 
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Figure 3: The graph showing schematically the asymptotic behavior of the n-th term of —E\s in the large-/?o 
approximation for n^> 1. 


Po approximation is known to lead to an O(Aqcd) uncertainty of the pole mass; on the other 
hand, the small size of the 0{a b s m) plus O(o.gm\ogas) corrections in the range fi ~ 10-30 GeV 
appears to be incompatible with the renormalon picture. 

Let us recall the estimate of the renormalon uncertainty in the large-/?o approximation (see 
e.g. PI )• Asymptotically the perturbative series of the IS 1 energy level, if expressed in the pole 
mass, behaves as 

E\s ~ const, x iias{n) x j x n \ f or n ^ i (4) 

It becomes minimal at n ~ n* = 2tt/(Pq as(/i)). The size of the term scarcely changes within 
the range n G (n* — y/ru, n* + y/rE); see Fig. |3]" We may consider the uncertainty of this 
asymptotic series as the sum of the terms within this range, since we are not sure where to 
truncate the series within this range: 


6E- 


is 



A, 


QCD- 


(5) 


The /r-dependence vanishes in this sum, and this leads to the claimed uncertainty. This ar¬ 
gument shows that when the relevant coupling constant as (/i) is small (corresponding scale /1 
is large), n* is large. Then each term of the series for n € (u* — y/n^, n* + y/n^) can become 
considerably smaller than Aqcd- 

According to this argument, the small size of the 0(a 5 s m ) plus 0(a 5 s m log 0 : 5 ) correction 
at certain scales does not generally lead to an uncertainty considerably smaller than Aqcd- 

1 Using the Stirling formula, one may easily find an approximate position n* of the minimum of the series. 

Then, by expanding around the minimum, one finds an approximate form n\n~ n « y/2i: n* exp[—n* + (n — 

n*) 2 * /(2?r*)] ~ y/2irnl exp(—n*) in the range |n — n*| < sjru. 
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While an error estimate should necessarily be more or less subjective, as long as the large- 
Pq approximation is valid, we should at least bear in mind how the theoretical uncertainty is 
estimated in this framework. Incidentally, based on the large-/£?o approximation, the MS mass 
extracted from the IS' energy level has an uncertainty of order Aq CD /( asm t ) 2 r\-/ AqCd/ opt rv 3 — 
10 MeV originating from the next-to-leading order renormalon contribution [2^. Thus, the 
above perturbative error of order 40 MeV is still significantly larger than this contribution. 

To conclude, we observe a much more stable prediction of the toponimn IS energy level 
when we use the MS mass instead of the pole mass. Considering this situation and the good 
agreement of the large-/3o approximation with the presently known corrections, we consider a 
theoretical uncertainty of the pole mass of order Aqcd~ 200-300 MeV to be legitimate. On 
the other hand, based on the argument in IQ, it is likely that the top quark MS mass can be 
extracted with an accuracy of order 40 MeV, once the 4-loop relation between the pole and 
MS mass is calculated. This number may be compared with the most recent estimate |21| of 
the experimental error (including some systematic errors) of 19 MeV in the determination of 
the top quark IS' mass, corresponding to a 3-parameter fit with an integrated luminosity of 
300 fb" 1 . 
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Appendix 

In this appendix we list the formulas we use to convert the energy level of the quarkonium IS 1 
state from the pole mass scheme to the MS mass scheme using the e-expansion H- 
The energy of the quarkonium IS state is given by 


M 1S = 2m po ie + E ls (m po i e , a s (n )) (6) 

as a function of m po i e and «s(/i) = a^ f (/r) in the pole mass scheme, where rif is the number of 
light quark flavors (rif = 4, 5 for the bottomonium and toponimn, respectively). Mass relation 
between the pole and MS masses is given by 


ttt po le \ 1 T d{ 


eas(m) 


7r 


+ d\ 


' £a>s(m) 


7r 


+ d 2 


' £as(m) 


+ d 3 


' £Q>s{rn) 


7T 


7T 


+ <?(£ 5 ) ,(7) 


where £ — 1 is the expansion parameter in the e-expansion, m = ^Ms( m Ms)’ — 4/3. 
The coefficients d\ and d 2 are obtained from the 2-loop El and 3-loop |9| mass relations **, 

**The same relation was obtained numerically before in \‘24\ in a certain approximation. 
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respectively, by rewriting them in terms of the coupling of the theory with rif flavors only. 
These are given by 


j _ 307 , 7T 2 , 7T 2 log 2 C 3 , 

32 3 9 6 1 

~ 13.4434- 1.04137/1/, 


71 


144 


7T 

18 


( 8 ) 


d 2 = 


8462917 652841 vr 2 695 vr 4 575 vr 2 log 2 


93312 


38880 


7776 


162 


22 vr 2 log 2 2 55 log 4 2 220Li 4 (|) 58 Cs 1439 tt 2 C 3 1975 Cs 


81 


162 


27 


+ 


27 


432 


+ 


216 


( 231847 991 tt 2 61 tt 4 11 tt 2 log 2 2 7r 2 log 2 2 log 4 2 

+Uf I ^ 7 ^-7777“ + —T-77- 1 -77-h 


23328 


648 


1944 


81 


81 


81 


iLLfi' 


27 


241 

72 


+ n-. 


{ 2353 13t r 2 , 7Cs 

V23328 + "~324 


+ 


54 


~ 190.391 - 26.6551 n f + 0.652691 nj, 


(9) 


with ( 3 = 1.20206 ■ • •, Li 4 (|) = 0.517479 ■ ■ ■. The third coefficient d 3 is not known exactly yet. 
In this paper we use its value in the large-/?o approximation [7j: 


, „ „ . & (42979 89tt 2 71 tt 4 317Cs 

d 3 ( large-/ 3 0 ) = ^ + -77777- + S3 


64 V 5184 ' 18 ' 120 

~ 3046.29 - 553.872 n f + 33.568 n f 2 - 0.678141 n f 3 . 


12 

2 


( 10 ) 


To achieve the renormalon cancellation between 2m po \ e and E x s{rn po \ e , as{p)) order by order 
in the e-expansion, we must use the same coupling constant «s(/i) in the series expansions of 
2m po i e and Ei$. Therefore, «s(m) is re-expressed in terms of a 5 (/i) as 

as(m) = as(fi)ll + (3 ° l0g ^ + ( ,]l l0g (^) + f3 ° lQg (iM/bp 


TT ) 

fo lQ g te) 5 A) A lo S 2 (i) ^0 Ip g 3 (m) \ / easil y 

32 32 8 I { tt I 


TT I 


+ 0(e 4 ) 


( 11 ) 


using the coefficients of the QCD /3-function: 


A> = 11 


2 n f 


Pi = 102 - 


38 71 / 


,, 2857 5033 nt 325 n f 2 

P2 = — -777 ~ + 1 


3 3 2 18 54 

Using Eqs. CD and (ED, we obtain the e-expansion for m po i e in terms of as (//), 

/ 4 


"ipoie = m x 1 + d n -i(lfi) £ r 


n =1 


^e(/i) y 


+ 0(s 5 ) , 


( 12 ) 


(13) 
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where the coefficients d n {l are functions of l, L = log(/i/m) which enter via Eq. (TTH) . 
The binding energy E ls (m po i e , a s (fji)) is given by 


As = -A(A 2 ™Poief> n+1 YdA) p„(£ ), (14) 

9 n=0 \ n J 

where L M = log [n/(C F as(n)m po i e )], and P„(L M ) are given by 

Po(L^) = 1, Pi(L^) — Po Lfj-\- Ci, 

P'iiL/j) = -0 q L^ -/?o + + 2^° Cl ) + ° 2 ' 

P 3 (L„) = i(3 0 3 L„ 3 + (-gft 3 + + |$ Cl ) 

+ [^1^0 — ^AA + IgA — ^/^0 C l + gAci + 2/3 0 C2^ A. + c 3 > (15) 


with 


Cl 

c 2 


97 11 

IT ~~ ~9~ 

1793 2917tt 2 


-n /; 


12 


216 



The terms which contain L M are determined by renormalization-group equation from lower order 
constants 2 - The Ci )2 are taken from [22J, C 3 is given in Eq. ©• To obtain the e-expansion in 
the MS scheme, we re-express the pole mass in Eus^poie, a s(l^)) by m and «s(/i) employing 
the mass relation Eq. m , which gives 


PlS ~ —g ^(h) 2 


m y £"+i (7AA p„(r„, y. +o(A), 


(17) 


n=0 


with = log [/r/(C'i?a 5 (/i)m)]. Using the e-expansions Eqs. (Td and m Mis is rewritten 


as 


Mis = 2mil + J2 d n -i(l^) e r 


= 2 m 


n= 1 
4 


(«M)-) _! 


a s (^) 2 m 


_n+1 


1 + ^ ^n-l(^) — 


\n=0 

27ras(/i) ~ , A / e 


'«A/YY 


Pn(Pfii 1 /j 


n =1 


9 


1 Pn —1 {Phi Y) ) ^ 


«Ah)Y 


(18) 


Setting the expansion parameter e = 1 in the final expression, the n-th order correction to Mi 5 
in the MS scheme is given by 2m x (as(/i)/7r) ?t d n _i(A) — (27ras(/r)/9) P ra -i(T M , Z M ) . 
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